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Abstract - The frequencies of a cryogenic 
sapphire oscillator and a hydrogen maser 
are compared to set new constraints on a 
possible violation of Lorentz invariance. We 
determine the variation of the oscillator 
frequency as a function of its orientation 
(Michelson-Morley test) and of its velocity 
(Kennedy-Thorndike test) with respect to a 
preferred frame candidate. We constrain the 
corresponding parameters of the Mansouri and 
Sexl test theory to 5 - f3 + 1/2 < 3.4 x 1CT 9 and 
(3 — a — 1 < 4.1 x which is of the same order 
as the best previous result for the former and 
represents a 50 fold improvement for the latter. 
These results correspond to an improvement 
of our previously published limits [Wolf P. et 
al., Phys. Rev. Lett. 90, 6, 060402, (2003)] by 
about a factor 2. We describe the changes of 
the experiment, and show the new data that 
lead to that improvement. 

The Einstein equivalence principle (EEP) is at the 
heart of special and general relativity jlj] and a corner- 
stone of modern physics. One of the constituent ele- 
ments of EEP is Local Lorentz invariance (LLI) which, 
loosely stated, postulates that the outcome of any local 
test experiment is independent of the velocity of the 
freely falling apparatus (the fundamental hypothesis 
of special relativity). The central importance of this 
postulate in modern physics has motivated tremen- 
dous work to experimentally test LLI pj. Addition- 
ally, nearly all unification theories (in particular string 
theory) violate the EEP at some level [2| which further 
motivates experimental searches for such violations of 
the universality of free fall Q and of Lorentz invariance 



The vast majority of modern experiments that test 
LLI rely essentially on the stability of atomic clocks 
and macroscopic resonators [E IS IE EiJ > therefore 



improvements in oscillator technology have gone hand 
in hand with improved tests of LLI. Our experiment is 
no exception, the large improvement being a direct re- 
sult of the excellent stability of our cryogenic sapphire 
oscillator. Additionally its operation at a microwave 
frequency allows a direct comparison to a hydrogen 
maser which provides a highly stable and reliable ref- 
erence frequency. 

Numerous test theories that allow the modeling and 
interpretation of experiments that test LLI have been 
developed. Kinematical frameworks [ill Il2| postulate 
a simple parametrisation of the Lorentz transforma- 
tions with experiments setting limits on the deviation 
of those parameters from their special relativistic val- 
ues. A more fundamental approach is offered by the- 
ories that parametrise the coupling between gravita- 
tional and non-gravitational fields (THe/x 0, Il3l IbH 
or XS fl5| formalisms) which allow the comparison of 
experiments that test different aspects of the EEP. Fi- 
nally, formalisms based on string theory IS HI have 
the advantage of being well motivated by theories of 
physics that are at present the only candidates for a 
unification of gravity and the other fundamental forces 
of nature. 

Owing to their simplicity the kinematical frame- 
works of [ill flij l have been widely used to model 
and inter pret many previous experiments testing LLI 
[E H 0. lid - ■ In order to compare our results 
to those experiments we will follow this route in the 
present work (an analysis of our experiment in the 
light of other test theories being relegated to a future 
publication). Those frameworks postulate generalized 
transformations between a preferred frame candidate 
E(T, X) and a moving frame S(t, x) where it is assumed 
that in both frames coordinates are realized by iden- 
tical standards (e.g. hydrogen masers for the time co- 
ordinates and sapphire rods for the length coordinates 
in our case). We start from the transformations of 12] 
(in differential form) for the case where the velocity of 



S as measured in E is along the positive X-axis, and 
assuming Einstein synchronization in S (we will be con- 
cerned with signal travel times around closed loops so 
the choice of synchronization convention can play no 
role): 
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with c the velocity of light in vacuum in E. Using 
the usual expansion of the three parametrs (a rs 1 + 
av 2 /c 2 + 0(c- 4 );6 ps 1 + f3v 2 /c 2 + £)(c- 4 );d « 1 + 
6v 2 /c 2 + 0(c- 4 )), setting c 2 dT 2 = dX 2 + dY 2 + dZ 2 
in E, and transforming according to (1) we find the 
coordinate travel time of a light signal in S: 
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where dl = yj dx 2 + dy 2 + dz 2 and 9 is the angle be- 
tween the direction of light propagation and the veloc- 
ity v of S in E. In special relativity a = — 1/2; P — 
1/2; 5 = and (1) reduces to the usual Lorentz trans- 
formations. Generally, the best candidate for E is 
taken to be the frame of the cosmic microwave back- 
ground (CMB) [H [13 with the velocity of the solar 
system in that frame taken as Vq w 377 km/s, deck 
w -6.4°, RA w 11.2h. 

Michelson-Morley type experiments 0, |i(J deter- 
mine the coefficient Pmm = (1/2 — p + 5) of the direc- 
tion dependent term. For many years the most strin- 
gent limit on that parameter was \Pmm\ < 5 x 10 -9 
determined over 23 years ago in an outstanding exper- 
iment 0| . Our experiment confirms that result with 
roughly equivalent uncertainty (3.4 x 10~ 9 ). Recently 
an improvement to |Pmm| < 1-3 x 10 -9 has been re- 
ported [2lJ . Kennedy-Thorndike experiments 0, H, 
measure the coefficient Pkt — (P—a—1) of the velocity 
dependent term. The most stringent limit [t| on |Prtt| 



has been recently improved from by a factor 3 to 
\Pkt\ < 2.1 x 10~ 5 . We improve this result by a factor 
of 50 to I Pkt I < 4.1 x 10~ 7 . Finally clock comparison 
and Doppler experiments [Til Il7l Y2fl\ measure a, cur- 
rently limiting it to |a + 1/2| < 8 x 10~ 7 . The three 
types of experiments taken together then completely 
characterize any deviation from Lorentz invariance in 
this particular test theory. 

Our cryogenic oscillator consists of a sapphire crystal 
of cylindrical shape operating in a whispering gal lery 
mode (see fig. 1 for a schematic drawing and [23L l2-l| 
for a detailed description). 
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FIG. 1: Typical relative frequency stability of the CSO - 
H-maser difference after removal of a linear frequency drift. 
The inset is a schematic drawing of the cylindrical sapphire 
oscillator with the Poynting vector P in the whispering 
gallery (WG) mode, the velocity v(t) of the cylinder with 
respect to the CMB, and the relevant angles for a photon 
in the WG mode. 

Its coordinate frequency can be expressed by v — 
m/t c where t c is the coordinate travel time of a light 
signal around the circumference of the cylinder (of ra- 
dius r) and m is a constant. From (2) the relative fre- 
quency difference between the sapphire oscillator and 
th e hy drogen maser (which realizes coordinate time in 
S 0) is 
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where i/q = m/(2Trr/c), v(t) is the (time dependent) 
speed of the lab in E, and tp is the azimuthal angle of 
the light signal in the plane of the cylinder. The peri- 
odic time dependence of v and 9 due to the rotation and 
orbital motion of the Earth with respect to the CMB 
frame allow us to set limits on the two parameters in 
(3) by adjusting the periodic terms of appropriate fre- 
quency and phase (see [13 f° r calculations of similar 
effects for several types of oscillator modes). Given the 
limited durations of our data sets (< 15 days) the dom- 
inant periodic terms arise from the Earth's rotation, so 
retaining only those we have v(t) = u + f! x R with 
u the velocity of the solar system with respect to the 
CMB, ft the angular velocity of the Earth, and R the 
geocentric position of the lab. We then find after some 
calculation 

Av/vq = P KT (Hsm\) + P M m (AcosA + £cos(2A) 
+CsinA + DsinAcosA + £'sinAcos(2A)) 

(4) 

where A = Qt+(j), and A-E and <j> are constants depend- 
ing on the latitude and longitude of the lab (« 48.7°N 
and 2.33°E for Paris). Numerically H w -2.6 x 1CT 9 , 
A « -8.8 x 1(T 8 , B « 1.8 x 1CT 7 , C-E of order 1(T 9 . 
We note that in (4) the dominant time variations of 
the two combinations of parameters are in quadrature 
and at twice the frequency which indicates that they 
should decorelate well in the data analysis allowing a 
simultaneous determination of the two (as confirmed 
by the correlation coefficients below). Adjusting this 
simplified model to our data we obtain results that dif- 
fer by less than 10% from the results presented below 
that were obtained using the complete model ((3) in- 
cluding the orbital motion of the Earth) . 

The cryogenic sapphire oscillator (CSO) is compared 
to a commercial (Datum Inc.) active hydrogen maser 
whose frequency is also regularly compared to caesium 
and rubidium atomic fountain clocks in the labora- 
tory 28). The CSO resonant frequency at 11.932 GHz 
is compared to the 100 MHz output of the hydrogen 
maser. The maser signal is multiplied up to 12 GHz 
of which the CSO signal is subtracted. The remaining 

67 MHz signal is mixed to a synthesizer signal at 
the same frequency and the low frequency beat at w 
64 Hz is counted, giving access to the frequency differ- 
ence between the maser and the CSO. The instability 
of the comparison chain has been measured and does 
not exceed a few parts in 10 16 . The typical stability 



of the measured CSO - maser frequency after removal 
of a linear frequency drift is shown in fig. 1. Since 
January 2003 we have implemented an active temper- 
ature control of the CSO room and changed some of the 
electronics. As a result the diurnal and semi-diurnal 
temperature variations during measurement runs (» 2 
weeks) were typically « 0.1° C in amplitude with best 
values < 0.04°, and longer and more reliable data sets 
became available. 

Our experimental data consists of two blocks, before 
and after the Jan. 2003 upgrade. The first block con- 
sists of eight sets of measured values of Av/vq taken 
between Nov. 2001 and Oct. 2002. Our previously 
published results |l£j are based on the first seven of 
these (3 to 10 days individual lengths, 37 days in total). 
The second includes 5 data sets taken between Jan. 
2003 and Apr. 2003 with significantly longer lengths 
(6 to 15 days, 56 days in total) and better tempera- 
ture stability. The sampling times To are generally 100 
s except for the first two sets for which To = 12 s and 
5000 s respectively. To analyze our data we simulta- 
neously adjust an offset and a rate (natural frequency 
drift, typically w 2 x 10~ 18 s _1 ) per data set and the 
two parameters of the model (3). In the model (3) we 
take into account the rotation of the Earth and the 
Earth's orbital motion, the latter contributing little as 
any constant or linear terms over the durations of the 
individual data sets are absorbed by the adjusted off- 
sets and rates. 

When carrying out an ordinary least squares (OLS) 
adjustment we note that the residuals have a signifi- 
cantly non-white behavior as one would expect from 
the slope of the Allan deviation of fig. 1 . The power 
spectral density (PSD) of the residuals when fitted 
with a power law model of the form S y (f) = kf^ 1 yields 
typically fi « — 1.5 to — 2. In the presence of non- white 
noise OLS is not the optimal regression method [2^.l30| 
as it can lead to significant underestimation of the pa- 
rameter uncertainties 29). 

An alternative method is weighted least squares 
(WLS) [13 

which allows one to account for non- 
random noise processes in the original data by pre- 
multiplying both sides of the design equation (our 
equation (3) plus the offsets and rates) by a weight- 
ing matrix containing off diagonal elements. To deter- 
mine these off diagonal terms we first carry out OLS 
and adjust the S y (f) = kf^ power law model to the 
PSD of the post-fit residuals determining a value of 
\x for each data set. We then use these [i values to 



construct a weighting matrix following the method of 
fractional differencing described, for example, in |29| . 
Figure 2. shows the resulting values of the two param- 
eters (Pkt and Pmm) for each individual data set. 
As previously reported a global WLS fit to the 
first 7 data sets yields \P M m\ = (1.5 ± 3.1) x 1CT 9 
and \P K t\ = (-3.1 ± 3.7) x 1(T 7 (la uncertainties), 
with the correlation coefficient between the two param- 
eters less than 0.01 and all other correlation coefficients 
< 0.07. Because of the qualitative difference of the 
data we decided to treat the recent data (since Jan. 
2003) independently, the global WLS adjustment for 
those 5 data sets yielding \Pmm \ = (5.5 ± 2.5) x 10 -9 
and \Pkt\ = (1.2±3.4) x 10~ 7 , again with insignificant 
correlation coefficients. 
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FIG. 2: Values of the two parameters (Pkt and Pmm) 
from a fit to each individual data set (blue diamonds), to 
the first seven data sets (green triangle, as previously pub- 
lished |lQ|), and to the last five data sets (red square, data 
since upgrade). The error bars indicate the combined un- 
certainties from statistics and systematic effects (see text 
for details). 

Systematic effects at diurnal or semi-diurnal fre- 
quencies with the appropriate phase could mask a pu- 
tative sidereal signal. The statistical uncertainties of 
Pmm and Pkt obtained from the WLS fit above corre- 
spond to sidereal and semi-sidereal terms (from (4)) of 
« 1 x 10~ 15 and as 5 x 10 -16 respectively so any system- 
atic effects exceeding these limits need to be taken into 



account in the final uncertainty. We expect the main 
contributions to such effects to arise from temperature, 
pressure and magnetic field variations that would af- 
fect the hydrogen maser, the CSO and the associated 
electronics, and from tilt variations of the CSO which 
are known to affect its frequency. Measurements of 
the tilt variations of the CSO show amplitudes of 4.6 
/irad and 1.6 /irad at diurnal and semi-diurnal frequen- 
cies. To estimate the tilt sensitivity we have intention- 
ally tilted the oscillator by « 5 mrad off its average 
position which led to relative frequency variations of 
f» 3 x 10~ 13 from which we deduce a tilt sensitivity 
of « 6 x 10~ 17 /xrad -1 . This value corresponds to a 
worst case scenario as we expect a quadratic rather 
than linear frequency variation for small tilts around 
the vertical. Even with this pessimistic estimate diur- 
nal and semi-diurnal frequency variations due to tilt 
do not exceed 3 x 10~ 16 and 1 x 10~ 16 respectively and 
are therefore negligible with respect to the statistical 
uncertainties. 

Before January 2003 temperature measurements of 
the CSO lab taken during some of the experimental 
runs show room temperature variations with ampli- 
tudes of 0.3 °C and 0.1 °C for the diurnal and semi- 
diurnal components which are reduced to 0.04 °C and 
0.01 °C on the actively temperature controlled elec- 
tronics panel. In december 2002 we implemented an 
active temperature stabilisation inside an isolated vol- 
ume (fa 15m 3 ) that included the CSO and all the as- 
sociated electronics. The temperature was measured 
continously in two fixed locations (behind the elec- 
tronics rack and on top of the dewar) the latter be- 
ing the reference point for the stabilisation servo loop. 
Using this system the diurnal/semi-diurnal tempera- 
ture variations observed were 0.1/0.08 °C in ampli- 
tude, for the worst data set and the measurement be- 
hind the electronics rack (on the dewar the fluctua- 
tions were always significantly smaller). The quoted 
numbers are the worst performance observed over all 
available data sets since Jan. 2003, in best cases they 
remained as low as 0.04/0.04 °C. The hydrogen maser 
is kept in a dedicated, environmentally controlled clock 
room. Measurements of magnetic field, temperature 
and atmospheric pressure in that room and the maser 
sensitivities as specified by the manufacturer allow us 
to exclude any systematic effects on the maser fre- 
quency that would exceed the statistical uncertainties 
above. When heating and cooling the CSO lab by 
k 3°C we see frequency variations of « 5 x 10~ 15 per 



°C. This is also confirmed when we induce a large si- 
nusoidal temperature variation (« 1.5°C amplitude). 
From the temperature measurements during the ex- 
perimental runs before Jan. 2003 we therefore deduce 
a total diurnal/semi-diurnal effect of w 1.5 x 10~ 15 
and « 0.5 x 10~ 15 respectively. For the more re- 
cent data (taking the worst temperature variations ob- 
served) we obtain w 0.5/0.4 x 10 -15 . We assume the 
pessimistic scenario where all of the diurnal and semi- 
diurnal effect is present at the neighboring sidereal and 
semi-sidereal frequencies which leads (from (3)) to un- 
certainties from systematic effects of ±5.8 x 10~ 7 on 
Pkt and ±2.8 x 10~ 9 on Pmm for the old data and 
±2.3 x 10~ 7 respectively ±2.2 x 10~ 9 for the data since 
Jan. 2003. We note that the phase of the perturbing 
systematic signal will vary over the course of our mea- 
surements due to natural causes (meteorology, daytime 
changes etc.) and to the sidereal/diurnal frequency 
difference so our final uncertainties are the quadratic 
sums of the above values and the statistical uncertain- 
ties from the WLS adjustment: Pmm < 3.4 x 10~ 9 
and \Prt\ < 4.1 x 10~ 7 for the new data. 

In summary, we have reported an experimental test 
of Lorentz invariance that simultaneously constrains 
two combinations of the three parameters of the Man- 
souri and Sexl test theory (previously measured indi- 
vidually by Michelson-Morley and Kennedy-Thorndikc 
experiments). Using our recent results (after Jan. 
2003) we obtain \S - /3 + 1/2| = 5.5(2. 5)(2.2) x 10~ 9 
which is of the same order as the best previous results 
HIHl, and \/3-a-l\ = 1.2(3.4) (2.3) x 10~ 7 which im- 
proves the best present limit |9( by a factor of 50 (the 
first bracket indicates the uncertainty from statistics 
the second from systematic effects). The measured de- 
viation of \5 — /3±1/2| from zero is significant at the la 
level (67 % confidence) and could indicate a deviation 
from Lorentz invariance. However, at the time being, 
we consider this to be most likeley a statistical artefact 
or an overlooked or insufficiently evaluated systematic 
effect. The experiment is running continuously and we 
expect that future data will help to clarify this situa- 
tion. As a result of our experiment the Lorentz trans- 
formations are confirmed in this particular test the- 
ory with an overall uncertainty of < 8 x 10 -7 limited 
now by the determination of a from Doppler and clock 
comparison experiments [lfil ll7T ] . This is likely to be 
improved in the coming years by experiments such as 
ACES (Atomic Clock Ensemble in Space H3) that will 
compare ground clocks to clocks on the international 



space station aiming at a 10 fold improvement on the 
determination of a. 
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